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Abstract. Quite detailed theoretical description of superfluid 3 He is possible on length scales that are much larger than the 
atomic scale. We discuss weak links between two bulk states of 3 He-B. The current through the weak link is determined 
by the bound states at the link. The bound state energies are spin split depending on the order parameters in the bulk. As a 
result, unusual current-phase relations with % state appear. For not too weak links, the order parameter in the bulk is modified 
because of the Josephson coupling. This leads to a stronger 7t state and to an additional current at constant pressure bias. The 
theoretical results are compared with experiments. 



INTRODUCTION 

In spite of more than thirty years of studies of superfluid 
phases of 3 He, there are several active research direc- 
tions. One of them is the study of Josephson phenomena 
in weak links between two volumes of bulk superfluid. 
Here we give a review of recent theoretical work on this 
subject QH^yM. We pay special attention to the role of 
bound quasiparticle states in the link, and explain how 
these determine the flow through the junction. In par- 
ticular we study the equilibrium current-phase relation, 
/((])), and the dc current at constant potential difference, 
hc(U). We compare the results with experiments. The 
experimental work has been reviewed in Ref. 1 5]. A the- 
oretical review with a somewhat different emphasis has 
been given in Ref. . 



WEAK LINKS 

Let us consider two volumes of superfluid that are con- 
nected by a weak link. Here "superfluid" means gener- 
ally a fermion superfluid, either superconducting metal 
or liquid 3 He. Our discussion of superconducting weak 
links is mainly aimed as an introduction for 3 He. There- 
fore we neglect all complications that arise from impurity 
or interfacial scattering in superconductors. Also we con- 
sider only superconductors of the conventional (s-wave) 
type. In the case of 3 He we limit to the superfluid B 
phase. 

The geometry of the weak link is depicted in Fig. ^ 
The figure shows one quasiparticle trajectory with mo- 
mentum p. In bulk superfluid such elementary excita- 
tions must have energy e larger than the energy gap A of 
the superfluid. In the weak link, however, there exists en- 




FIGURE 1. Sketch of a weak link between two bulk super- 
fluids. The two bulk superfluids have phases 0^ and (])#. One 
quasiparticle trajectory with momentum p is shown. When its 
energy is within the energy gap, the quasiparticle cannot escape 
to the bulk but is Andreev reflected as hole with essentially the 
same momentum. This in turn is Andreev reflected as particle 
on the other side of the link. Thus particle and hole type exci- 
tations bounce back and forth along the same trajectory. 

ergy eigenvalues within the energy gap, |e| < A. In such 
a state the quasiparticle cannot escape to the bulk but is 
Andreev reflected. In Andreev reflection a particle type 
excitation is converted to hole type and vice versa, with 
essentially unchanged momentum [7]. Thus the particles 
and holes travel the same trajectory back and forth, re- 
spectively. In one cycle, one Cooper pair, or equivalently 
two particles, are transmitted through the weak link. 

The bulk superfluid states are characterized by phases 
§l and §r- The energies of the bound states depend 
essentially on the phase difference (]) = §r — (j)^ |8|. For 
superconductors this is depicted in Fig. |2] The states 
appear as pairs with positive and negative energies, ±|e|. 
Understanding that dt/d§ plays a similar role as the 
group velocity, we can assign the states as propagating 
to the right (8 = +1) and to the left (8 = — 1). These two 
states are not independent, however. Based on what was 
stated about Andreev reflection above, if a state with 8 = 




FIGURE 2. Bound state energies £g(<|>) in a superconducting 
point contact. As (]) increases, the states with transport to the 
right (8 = +1) have increasing energy (solid lines), and the 
states with transport to the left (8 = — 1) have decreasing energy 
(dashed lines). The states at energies |e| > A form a continuum. 

+ 1 is occupied, the corresponding 8 = — 1 state must be 
empty. And vice versa. More generally, the occupations 
f± have to satisfy /+ + /_ = 1 . Really the states with 
positive and negative energies are the same bound state, 
whose energy relative to the ground state is the absolute 
value |e|. Anyhow, we find it very useful to consider a 
bound state as a superposition of positive and negative 
energy states. 

In the simplest case of a point contact, the bound 
state energies have the simple analytic expression 8 = 
± A cos | 0] . By a point contact, or a pinhole in case 
of 3 He, we mean a weak link whose all dimensions are 
small compared to the superfluid coherence length ^o- 

It is an interesting observation that to a large extent the 
current th roug h a weak link is determined by the bound 
states OH In the following we are interested 
in equilibrium or only small nonequilibrium, where the 
chemical potential difference ft-/)jE(/«A~ ksT c . 
In this case the current has the form 

«> 

Here / ; g is the occupation of the state. In equilibrium it 
reduces to the Fermi distribution = [exp^^/ksT) + 
1] . The index i indicates different channels, which 
correspond to different directions and locations of the 
bound state (Fig. 0. In a point contact of area A, there 
are M — k\A /An channels, where kf is the Fermi wave 
vector. In the normal state the current is / = MU j%h [12]. 

Kulik and Omel'yanchuk derived in 1975 the follow- 
ing formula for equilibrium current trough a supercon- 
ducting point contact [13] 

MA d) Acos$ 

7= sin-tanh (2) 

Ti 2 2k B T 

Their original derivation does not give much clue to un- 
derstand the result. Now this result can be straightfor- 
wardly understood based on the formulas given above: 
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FIGURE 3. Example of bound states energies e ; g no ((|)) in a 
3 He weak link. Compared to Figure|2| the bound state energies 
have smaller slopes and they are spin split (a = ±1). Figure 
reprinted from Ref. y[. 

all channels have the same bound state energies e = 
±Acos | and using the current formula Q and the Fermi 
distribution gives the result Q. 

As the temperature increases towards the superfluid 
transition temperature T c , the Kulik-Omel'yanchuk re- 
sult (|2j reduces to sinusoidal form 

7 = / c sin<t) (3) 

with critical current I c = MA 2 /AhksTc. 

3 HE WEAK LINKS 

Most of the discussion of superconducting weak links 
applies also to weak links in 3 He. There are a few dif- 
ferences. In 3 He the Cooper pairs are in p-wave states, 
instead of the s-wave state considered above. Any scat- 
tering breaks these pairs. Thus the superfluid state is al- 
ways suppressed near walls. This suppression also af- 
fects the bound states energies in a weak link, as de- 
picted in Fig. [3] Compared to the superconducting case, 
the slopes de/d§ are smaller in magnitude. Another dif- 
ference is that there is spin splitting of the energy states 
(o = ±1). Whereas the states in the superconducting 
case were doubly degenerate, this degeneracy is gener- 
ally broken in 3 He. 

Because of the spin splitting the current formula Q 
has to be generalized to the form 

^iLLLL^f^ (4) 

where a factor of 2 is replaced by spin summation. There 
is also additional summation over index n, which takes 
into account that several bound states can occur at a given 

The spin splitting can lead to current-phase relation- 
ships that differ essentially from the standard sinusoidal 
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FIGURE 4. Current-phase relationships calculated for a pin- 
hole. The left hand panel is for parallel n vectors on the two 
sides of the junction. The right hand panel is for antiparal- 
lel n's that are perpendicular to the wall. At high tempera- 
ture the curves are sinusoidal At lower temperature the 
curves become skew in the parallel case resembling the Kulik- 
Omel'yanchuk result 0. The antiparallel case has negative 
critical current I c , or equivalently, is shifted by a phase differ- 
ence jr. At very low temperature it develops an additional kink 
that is known as 71 state. Figure adapted from Ref. 1 2]. 



form OJl |14]. Consider, for example, the case where the 
spin split states are shifted relative to each other by a 
phase difference rj jc. It follows that the leading sinu- 
soidal components of the current from the two spin states 
[°c sin((j) + \|/ a )] cancel each other. What is then left are 
higher harmonics sin(n0 +\|/ ), n = 3,...]. Unusual 
current-phase relations are confirmed by calculations be- 
low. 

The Cooper pairs of superfluid 3 He have p-wave form 
fisll . There are three orthogonal p-wave states: p x , p y , 
and p z . The spin state is triplet. There are three triplet 
states, which conventionally are chosen as — ft + ll> 
' TT +' ll> an d — Tl + IT- Here we concentrate on the 
B-phase, where p x pairs have the first spin state, p y pairs 
the second and p z pairs the third. However, the spin 
coordinate axes (*',/, z') are rotated relative to the orbital 
axes (x,y,z). The rotation angle is fixed to 104°, but the 
axis n of this rotation can vary. 

In hydrodynamics of the B phase the mass ((])) and spin 
(n) degrees of freedom are independent. The weak link 
acts as a nonlinear element that couples mass and spin, 
as will be demonstrated below. 



Isotextural Theory 

Some current-phase relations calculated for a 3 He-B 
pinhole are shown in Fig.0] Two cases are shown. First, 
the spin rotation axes n on both sided of the junction are 
the same, and second, the rotation axes have opposite di- 
rections perpendicular to the wall of the pinhole. The for- 
mer case resembles the Kulik-Omel'yanchuk result 0: 
at high temperatures 7(0) is sinusoidal but at lower tem- 




FIGURE 5. Measured current-phase relationships for a 65 x 
65 array of 100 nm diameter apertures in 50 nm thick wall Hal . 
In cooling through T c , the system randomly freezes to either H 
or L state. Only the relative value of the phase is determined 
experimentally. Figure adapted from Ref. Ilfll . 



peratures it becomes skew towards = 71. The case with 
antiparallel n's is also sinusoidal at high temperatures, 
but has negative I c , (Such a case is called a tc junction in 
superconductivity.) At temperatures around 0.1 7^ it de- 
velops to a TC state. This means that the slope of 7((j)) is 
positive at both and TC. 

Prior theoretical calculations, the tc state had been ob- 
served experimentally. The experimental results of Ref. 
1161 are shown in Fig.|5] The experimental 7(0) is sinu- 
soidal at high temperature but develops a tc state at low 
temperature. Moreover, two different metastable states 
with different sets of 7(0, r) curves was observed. These 
are called H and L, corresponding to high and low criti- 
cal current. The two states appeared randomly in cooling 
through T c , but the state selected remained stable during 
the stay in the superfluid state. 

Let us consider a weak link under a constant chemi- 
cal potential difference U = pi — PR- According to the 
Josephson relation 



c/0 
dt 



2U 



(5) 



the phase increases linearly in time. This implies that 
the bound state energies (Fig. [3j are continuously shift- 
ing up or down. During shifting the occupations of these 
states can change in collisions with other quasiparticles. 
Since the collisions are rare, this thermalizes the occu- 
pations only at small U <n/x, where x is the scattering 
time. At larger bias, the occupations of the bound states 
remain essentially constant during the shift from —A to 
+A, or from +A to —A. The occupations of these states 
are then fixed to the thermal equilibrium occupations at 
the starting energies, —A or +A, respectively. 

It is now possible to calculate the current using Eq. 
@. The resulting time-averaged current 7d c is plotted in 
Fig-El At large bias (U ^> h/i) the occupation are deter- 
mined by gap edges and the current therefore becomes 
independent of the bias U . At small bias (U <C h/x) the 
scattering has time to preserve a nearly thermal distribu- 
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FIGURE 6. Average current vs. pressure bias. The data 
points are experimental results from Ref. fl7ll . and the lines 
are theoretical results from Ref. |4]. The pressure bias of 1 Pa 
corresponds to U = 4.7 x 1CP 3 kgT c . Figure reprinted from Ref. 
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FIGURE 7. Average current vs. pressure at low pressure 
bias. The data points are experimental results from Ref. 1 18] for 
the same 65 x 65 array of apertures as the equilibrium currents 
in Fig. [5] The lines drawn on the right hand side panel are the 
results of isotextural theory, and fitted with a single parameter 
(the relaxation time) as in Fig. |6| According to isotextural 
theory, the same lines should fit also the H state, which clearly 
is not the case. Figure adapted from Ref. Q]. 



tion, and the deviation from equilibrium distribution as 
well as the current is linear in U. Note that the equilib- 
rium current (Fig.|4ji is oscillating, and therefore does not 
contribute to the time-averaged current. 

Fig-IU shows also experimental data by Steinhauer et 
al. 11*711 . The theory can reasonably be fitted to this data 
using the scattering time as the only fitting parameter. 
The agreement is surprisingly good taking into account 
that the single aperture in the experiment (dimensions 
7.8 /jm x 0.27 fjm) is large compared to the superfluid 
coherence length cjo = 77 nm, whereas the pinhole theory 
makes just the opposite assumption. 

The current at finite bias was also measured for the 
same array of apertures as the equilibrium current in 
Fig. 13 The results are shown in Fig. [7] The current is 
different for the H and L states. The current in the L state 
is in reasonable agreement with theory using again the 
relaxation time as the only adjustable parameter. 

As a summary thus far, we can say that theory and ex- 
periment are in reasonable agreement. Both show bista- 



bility (H and L state), % states, and similar dissipative 
currents Id c (U). The temperature dependencies agree as 
well the order of magnitudes of all quantities. 

In spite of this success, closer examination reveals 
some problems in the theory. 1) Only one of theoreti- 
cal bistable states shows 7t state (Fig. 0}. More detailed 
analysis reveals that this problem cannot be removed by 
considering other configurations of n's [2]. 2) The the- 
oretical K state occurs at a much lower temperature and 
is much weaker than in experiments (compare Figs. @] 
and|5}- 3) The theory does not predict any difference in 
hc(U) between the H and L states. This is because dif- 
ferent n configurations mainly shift the bound states en- 
ergies e,g„ a along the § axis but do not much affect the 
shape of the e,§„ ((j)) curves. In running (j) this only af- 
fects the instantaneous current 7(f), but the average 7<j c is 
unchanged. 



Anisotextural Theory 

All the problems above can be explained by a single 
new concept. If the Josephson coupling is strong, it can 
change the spin-part of the order parameter on both sides. 
We call this anisotextural effect, since the n orientations 
are commonly known as texture. 

As an example consider the case of parallel n's. Then 
there is no spin splitting and the bound state at phase 
difference <|) = Jt lies at zero energy, e = 0. The thermal 
occupation of this doubly degenerate state is / = 1/2. If 
we now allow spin splitting, the e = state can split into 
two with positive and negative energies (Fig. [3}- Taking 
into account that the occupation of the negative energy 
state is larger than that of the positive energy state, this 
leads to a reduction of energy. Thus the spin splitting can 
take place spontaneously. 

The anisotextural effect as a function of phase is 
shown in Fig. [8] Starting from the completely spin- 
symmetric situation at = 0, the n texture changes spon- 
taneously as (j) increases towards K. The reduction of the 
Josephson coupling energy Fj is associated with a change 
in /((()) since 



Thus a 7t state develops if the reduction of Fj is sufficient 
to produces a local minimum of Fj((j)) at (]) = K. 

An additional contribution to energy arises from the 
bending of n (right hand panel of Fig. [§). In a simple 
model the total energy is written as 

F^^Fj^^ + ^Kjd^l^l 2 . (7) 

The first term is the Josephson coupling energy that 
depends on (j) and on the tilting angle T|o of n at the weak 
link. The second term is a gradient energy for the tilting 




FIGURE 8. Principle of the anisotextural Josephson effect. 
The panels on the left and right hand sides depict an array 
of apertures and the configurations of the spin-rotation axis n. 
Assuming there is no spin structure (n = constant) at the phase 
difference (j) = 0, the n texture changes spontaneously when 
<j) increases to Ji. This leads to a reduction of energy and to a 
positive slope of /((])) at <]) = n, as shown by the change from 
dashed lines to solid lines in the current and energy plots. 




FIGURE 9. Current-phase relationships calculated using the 
anisotextural model Q. The theory contains one free parameter 
T|oo, and the result should be compared with the experimental 
one in Fig. [5] Figure adapted from Ref. 0]. 



angle in the bulk. In comparison to experiments, the only 
free parameter is the tilting angle r|oo on one side far away 
from the junction, which is difficult to calculate because 
of the complicated shape of the experimental cell. Using 
the freedom to adjust r|oo, it is possible to generate the 
current-phase relations shown in Fig. [9] This should be 
compared with the experimental curves in Fig. [5] 

Next we consider the anisotextural effect at constant 
bias U. Finite U means that the phase changes at constant 
rate [Eq. Q], As a result of the anisotextural effect, 
the n texture oscillates at the angular frequency co = 
2U/h. This oscillation of the spin structure generates 
spin waves, which radiate out of the junction. 

In order to make an quantitative theory, we need to add 
one more contribution to the energy Q so that the total 
energy is 



F[r\,S]=Fj(i] , 



\K$d\\VT\\ 2 



(8) 



The last therm is the energy associated with net spin 
density S, where y is the gyromagnetic ratio and % the 



magnetic susceptibility 1 19]. Here S and r| are conjugate 
variables. Writing the Hamiltonian equations, this leads 
to a wave equation describing the spin waves and to a 
boundary condition for them. It turns out to be a quite 
standard radiation problem. The radiated power has the 
frequency dependence 



^rad ' 



or 



1 + (r Co/c) 



2 ' 



(9) 



where ro is the radius of the Josephson array, and c 
the spin wave velocity. This expresses that there is little 
radiation for wave lengths longer than the source, but 
more at shorter wave lengths. 

The dissipated power now has to come from the dc 
current, P = UI& C . Thus there is additional dc current due 
to the spin wave radiation. Including also the prefactors 
we get for it the expression 



Aic,rad — 



CO 



frad _ 2 [JsM°°)\ 

U h 4ncK l + (cor /c) 2 



(10) 



Here the equilibrium spin current 7 sp (r| o) plays a role 
of a coupling constant between mass and spin variables. 
Most interestingly, this coupling is different for the the- 
oretical H and L states. Putting in numbers we see that 
the theory explains approximately one third of the ob- 
served H-L difference shown in Fig. [5] (The comparison 
is shown in Ref. [3].) This is not completely satisfactory, 
but we must remember that all parameters of the theory 
were fixed before the comparison. 

Models to explain the observed (U ) have also been 
suggested in Ref. 11811 . These models differ essentially 
from the present ones. In particular, an A-phase-like 
distortion of the order parameter is associated with the 
linear part and a quasiparticle radiation mechanism with 
the nonlinear part of I& Q (U) . 



CONCLUSION 

We have seen that coupling between mass and spin vari- 
ables is essential for understanding the Josephson phe- 
nomena in superfluid 3 He. With spin splitting of bound 
states together with the anisotextural effect it is possi- 
ble to explain most of the experimental results of both 
/((])) and Idc(U). The theory presented above uses the 
pinhole approximation. For more accurate results it is 
necessary to consider larger apertures. This would re- 
quire selfconsistent determination of the order parameter 
within the aperture, which is calculationally demanding 
and has been done only in limiting cases l20H2lll . We ex- 
pect that this might lead to qualitatively same type of re- 
sults as the anisotextural effect for pinhole arrays, in par- 
ticular, to enhanced 7t states and to a texture dependence 



of Idc(U). This might explain the remaining differences 
between theory and experiment, and the experimental re- 
sults in single apertures 1 22[ . 

The anisotextural phenomena depend strongly on 
many parameters such as array size, geometry of the ex- 
perimental cell and magnetic field. We hope that these 
could be tested in future experiments. 
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